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ABSTRACT 

Context. Circumstellar debris disks are important because of their connection and interaction with planetary systems. An efficient 
way to identify these systems is through their infrared excess. Most studies so far concentrated on early-type or solar-type stars, but 
less effort has gone into investigating M-dwarfs, which pose the additional problem that their mid-infrared colors are poorly known. 
Aims. We characterize the infrared photometric behavior of M-dwarfs and search for infrared excess in nearby M-dwarfs taken from 
the volume-limited RECONS sample, concentrating on mid-infrared wavelengths corresponding to warm (> 100 K) dust. We then 
check whether the population of these late-type stars has a significantly different fraction of infrared excess compared to earlier-type 
stars. 

Methods. We used data recently released from the WISE satellite, which provides the most sensitive mid-infrared all-sky survey to 
date. Our sample consists of 85 sources encompassing 103 M-dwarfs. We compared this sample to Spitzer data and matched it to the 
2MASS catalog. We derived empirical infrared colors from these data and discuss the intrinsic spread within these colors as well as 
the errors from WISE and 2MASS. Based on this, we checked the stars for infrared excess and discuss the minimum excess we would 
be able to detect. 

Results. Other than the M8.5 dwarf SCR 1845-6357 A, where the mid-infrared excess is produced by a known T6 companion, we 
detect no excesses in any of our sample stars. The limits we derive for the 22/jm excess are slightly higher than the usual detection 
limit of ~ 10-15% for Spitzer studies, but including the 12/jm band and the [12] - [22] color in our analysis allows us to derive 
tight constraints for the fractional dust luminosity L^ustlL-,,. Comparing our findings to earlier-type stars, we show that this result is 
consistent with M-dwarf excesses in the mid-inrared being as frequent as excesses around earlier-type stars. The low detection rate of 
01q q% we derive for our sample in that case is an age effect. We also present a tentative excess detection at 22yum around the known 
cold debris disk M-dwarf AU Mic, which is not part of our statistical sample. 

Conclusions. There is still no clear detection of a mid-infrared excess around any old (>30 Myr) main-sequence M-dwarf. It is 
unclear whether this is due to a different dust evolution around M-dwarfs compared to earlier-type stars or whether this is an age 
effect combined with the difficulties involved in searching M-dwarfs for infrared excesses. A significantly larger sample of well- 
studied M-dwarfs is required to solve this question. At the same time, their behavior at longer wavelengths, which are sensitive to 
colder dust, needs further investigation. 
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1. Introduction 



Dust around main-sequence stellar systems has become a very 
well-studied phenomenon since it was first seen. The first de- 
tection was made aroun d Vega, an AO s tar, th rough the detec- 
tion of infrared excess (lAumann et al.L |1984|). A lot of stars 
showing this "Vega phenonenon" were subsequently detected 
and imaged through the all-sky survey performed by IRA S at 
25, 60 and 100 urn (IZuckermani 120011; iRhee eTaLl l2007h. as 
well as with pointed observations using ISO (|Decin et al.V2003'; 
Spang ler et al., 2001) at 60 and 100 yum, S pitzer (Rebull et al., 
."2008t (Carpenter et all l2009t iKoerner et al.1 120 10) at 24 and 70 
lum, and Herschel (most recently by Acke et aTri2012.) . Recently, 
the all-sky AKARI satellite was able to m ake detections at 9 and 
18;um wavelength dFuiiwara et al.Ll2009h . 

It was quickly realized that the infrared excess must stem 
from dust that surrounds the star in a disk. It was further- 
more understood that micron-sized dust cannot persist in orbit 
around a star for very long due to effects such as radiation pres- 
sure, Poynting-Robertson drag, and stellar wind pressure. Thus, 
the dust has to be permanently regenerated in situ. The stan- 



dard interpretation of this is that the debris traces underlying 
planetesimal belts, where d ust is generated through collisions 
( Bac kman & Parescel[T993h . The same effect happens in the so- 
lar system's asteroid belt and is responsible for zodiacal light. 
The difference to the solar system is that these belts are mostly 
collisionally dominated, producing a much stronger infrared ex- 
cess signal than the solar system debris disk does. 

The first detection of this phenomenon was made around 
an A-star. Because these bright stars are easiest to observe, 
the detection of an infrared excess is relatively easy for A, F, 
G, and nearby K-stars. Late-type star and specifically ]VI-dwarf 
infrared excesses are significantly harder to detect. First, they 
are much fainter, and secondly, the Rayleigh-Jeans approxima- 
tion does not hold in the mid-infrared for M-dwarfs. Molecular 
bands dominate their spectra at visible and near-infrared wave- 
lengths. Because of this, the expected scatter in their mid- and 
far-infrared colors is larger compared to early-type stars. Most 
searches for in frared excess concentrate on early-type or solar- 
type s tars (e.g.. lSu et al 1 12006t iBrvden et al.L l2009t iMoor et all 
1201 Ih . 
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A sec ond effect is that the fre quency of debris disks declines 
with age ('Carpenter et al.L l2009l) . Vega has an age of only ~450 
Myr ( Yoon et al., 2010), and many studies for infrared excesses 
around early-type stars include very young targets (e.g. lSu et all 
|2006|). Very few young main-sequence M-dwarfs are known; one 
of them is the well-st udied debris disk system AU Mic (e.g., 
iFitzgerald et al.Ll2007h . This is because of the intrinsic faintness 
of M-dwarfs, which only allows detecting them if they are close 
to the solar system. M-dwarf debris disk surveys are therefore 
usually bound to target mostly old stars, where the excess fre- 
quency is expected to be lower Of course, many TTauri stars are 
M-dwarfs with strong mid-IR excess from primordial disk emis- 
sion and there is also growing evi dence that primordial d isks last 
longer around lower-mass stars dCarpenter et all l2006h . but no 
mid-infrared excesses attributable to debris have been found ex- 
cept around AT Mic, like A U Mic a member of the ~ 12 Myr 
p Pic toris moving group (Gaut ier et al.L l2007t iPlavchan et al.L 
|2009|) . Few detections of cold, Kuiper-belt-like debris disks have 
been made in the sub-mm re gime for young M-dwarfs (iLiu et al.L 
I2004t iLestrade et al.L |2009|) . at a seemingly lower rate than for 
earher-type stars. It has been pointed out that there is an ap- 
parent lack of debris disks around older M-dwarfs in general 
(IPlavchan etan.l2005h . 

On the other hand, M-dwarfs make up the majority of stars 
in our galaxy, which is a good reason to study them. Because 
of its magnitude limits and the intrinsic faintness of M-dwarfs, 
IRAS can only give access to the very closest M-dwarf sys- 
tems. The recently flown WISE satellite does not have the far- 
infrared capabilities of IRAS, but goes significantly deeper in the 
12- and ll-jjm bands, giving access to more M-dwarf systems 
(IWri ght et al., 2010). M-dwarf studies going this deep have been 
possible earlier only with dedicated pointed obs ervations, for ex- 
ample using the Spitzer space telescope (.Werner et al.Ll2004h . 

Theoretical studies of the behavior and evolution of dust 
around M-dwarfs suggest that the timescales and frequencies 
involved in M-dwarfs might be consi derably differe n t from 
what we observe in earlier- type stars (iPlavchan et al.L l2005h . 
Observationally, there is a lack of solid evidence for a differ- 
ence in either the frequency or the timescales to disperse debris 
disks when comparing M-dwarfs to earlier-type stars. No stud- 
ies were able to prove a statistically sign ificant difference be - 
tween A-stars and solar-type stars, either (iTrillin^ et al.L l2008h . 
The question whether the processes involved in the evolution of 
debris disks work similar over a wide range of stellar luminosi- 
ties or whether there are significant differences is so far unre- 
solved. 

We describe the data we acquired from WISE and 2MASS 
in section 2. In section 3, we describe the analysis of the data, 
our efforts to detect mid-infrared excesses and the limits we can 
derive from this on the excess at different wavelengths and on the 
dust surrounding the M-dwarfs in our sample. In section 4, we 
compare our study to other results from other M-dwarf surveys 
as well as surveys of earlier-type stars. We discuss implications 
in section 5 and draw our conclusions in section 6. 

2. Source data from RECONS, WISE, and 2MASS 

The WISE satellite was launched on December 14, 2009 and be- 
gan its operations in January 2010, scanning the complete sky in 
four waveba nds centered at 3.4, 4.6, 12, and 22 fim, named wl 
through w4 (Wright et al.','2010'). On April 14, 20 11, a prelimi- 
nary data release covered about 57 percent of the sky. On March 
14, 2012, the all-sky catalog was released. The WISE satellite 
is significantly more sensitive than the IRAS satellite was, and 



is also more sen s itive t han the more recent AKARI sky survey 
(iMurakami et al.L l200'7h . Its 5-cr detection limits are 16.4, 15.4, 
11.1, and 7.8 magnitudes, respectively. The WISE beams have a 
FWHM of 6.1", 6.4", 6.5", and 12", respectively. This means 
that close binary systems cannot be resolved, an aspect we will 
have to take into account. 

The 2MASS survey covered the entire sky in the J, H, and 
Ks wavebands d uring the years 1997 th rough 2001 using ground- 
based facilities dSkrutskie et aUl2006h . Its sensitivity as well as 
its resolution are better than the WISE limits. The errors in the 
2MASS magnitudes stem mainly from atmospheric effects and 
cahbration uncertainties. 

2.1. Stellar sample 

Because of the intrinsic faintness of late-type dwarfs, and be- 
cause of the magnitude limit of WISE, which is still relatively 
strict at 7.8 magnitudes in the longest waveband, a study of the 
mid-infrared behavior and excesses has to focus on nearby M- 
dwarfs. Very late M-dwarfs (M7, M8, M9) are too faint for the 
WISE 22 jjm band even at distances of only 10 pc, while early- 
type M-dwarfs (MO, Ml) can be usefully tracked out to distances 
of 25 pc or more. Because our aim is to have a sample not bi- 
ased heavily toward early-type M-dwarfs, and because M8 and 
M9 dwarfs are intrinsically rarer than MO and Ml dwarfs (the 
IMF peaks around M3-M4, Bastian et al.Ll2010l) . we concentrate 
on the immediate solar neighborhood, omitting many early-type 
M-dwarfs that would be within reach of WISE, but are beyond 
lOpc. 

One of the best and most complete samples of the immediate 
solar vicinity is the list of the 100 nearest stellar systems by the 
RECONS team dJao et al., 2005; Henry et si., 2006) . Including 
eight systems that were pushed out of the top 100 list, it con- 
sists of 108 systems encompassing 139 stars (not counting plan- 
ets and white or brown dwarfs, which also appear in the cata- 
log). The vast majority of these stars are M-dwarfs (110). One 
of the advantages of the RECONS sample is that in addition to 
giving accurate positions and proper motions (required because 
these very near objects move significantly on sky between the 
2MASS and WISE surveys), it also gives accurate parallaxes, 
which can be used to determine absolute luminosities. V mag- 
nitudes are also provided. RECONS can thus give a clean, very 
complete volume-limited sample, which is ideal for our purpose. 
It was also readily available at the time we started our analy- 
sis, in contrast to, for example, the Reid 8pc sample (Reid et afl, 
2007), which would have been another option. The number of 
M-dwarfs in the Reid sample is slightly larger. However, the sta- 
tistical analysis we aimed to perform can be performed equally 
well with the RECONS sample, and the slightly larger number 
of M-dwarfs in the Reid sample would not change our discussion 
significantly. 

2.2. Observational data 

All M-dwarfs in the RECONS sample can be seen at the ex- 
pected positions in the WISE images. Because some of the stars 
are confused with nearby earlier stars and/or tight binary M- 
dwarf companions, there are however only 90 individual sources 
for the 110 M-dwarfs in the catalog. Out of these, five sources 
containing seven stars (GJ 783 B, GJ 338 A, LP 771-095 A/B/C, 
GJ 725 A and GJ 725 B) had to be removed from our analysis 
because of poor data quality, mostly from confusion with nearby 
sources. One more star (GJ 752 B) had to be excluded from the 
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Table 1. Excluded stars 



Name 



reason for exclusion 



GJ 783 B 
GJ 338 A 
LP 771-095 A/B/C 



GJ 725 A/B 



GJ 752 B 



confusion witfi nearby (7") K2.5 primary 
confusion witii nearby (17") K7 secondary 
WISE catalog shows multiple sources at ex- 
pected position, but none of them matches 
the star 

close binary system (14") that has very bad 
S/N ratio due to mutual contamination 
excluded from analysis of w4 band, not de- 
tected in WISE w4 



Notes. Stars that were excluded from the analysis for reasons of data 
quality. 
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Fig. 1. WISE [3.4]-[4.6] color (red stars), [3.4]-[12] color (black 
diamonds) and [3. 4] -[22] color (blue triangles), plotted against 
the 3 .4 yum absolute magnitude. Both the [3 .4] - [ 1 2] and the [3 .4] - 
[22] color show the expected trend (they become redder for late- 
type stars) and can be fitted with a polynomial. The [3.4]-[4.6] 
color is expected to show a similar trend and furthermore be 
smaller than the other two colors for any spectral type. This 
is clearly not the case for the bright end of our sample (above 
magnitudes of ~6 or absolute magnitudes of ~8-9). We there- 
fore omit the second channel in our analysis. 



analysis of the 22 fim band because it could not be detected at 
this wavelength. A summary of the excluded sources and the 
reasons for exclusion is presented in Table [T] This leaves 103 
M-dwarfs in 85 individual WISE sources (70 single, 12 double 
and 3 triple systems) for our analysis (one single star less for the 
w4 band). 

Table |2] shows the observational data for these 85 sources. 
The WISE data is complemented with magnitudes from the 
2MASS survey and spectral types from the literature as well as 
V band data and distance information from RECONS. If multi- 
ple stars were present in the WISE beam, the V magnitudes were 
calculated by combining the single-source V band magnitudes 
from RECONS. In the single case that stars were unresolved 
in WISE, but resolved with 2MASS (GJ 896 A/B), the 2MASS 
magnitudes were calculated in the same way. In all other cases, 
the sources were confused in 2MASS as well. 




7 6 5 

SPITZER 24 urn magnitude 

Fig. 2. Compariso n of WISE and Spit zer data for the sources 
in common with Gautier et al. 1 (l2007h . The red line shows a 
weighted least-squres fit for the offset between the Spitzer and 
WISE magnitudes, which is 0.083 mag with the WISE magni- 
tudes being fainter This offset stems from the different wave- 
length passbands (centered at 22 //m vs. 24 //m). The WISE data 
match the Spitzer data very well with only a few outliers at the 
faint end. 



For the purposes of this paper, the data from three of the four 
WISE bands were used: The first, third, and fourth. The second 
band, centered at 4.6 yum, was omitted. The reason is that while 
the w4 band is just about able to pick up the faintest and far- 
thest dwarfs in our sample at magnitudes beyond 8, the WISE 
measurements for the earlier bands run into saturation at magni- 
tudes brighter than ~6. While the first WISE band behaves quite 
well and high-quality magnitudes can be recovered, the magni- 
tudes of the second band are prone to significant errors for bright 
stars. Because of this, we did not trust the 4.6yum measurements 
for magnitudes smaller than ~6 (in our case, absolute magni- 
tudes smaller than ~8-9), which means large parts of our stellar 
sample. We illustrate this in Figure [T] Because of this, we fo- 
cused on the 3.4, 12, and 22 jum bands. 

2.2.1. WISE data quality 

The WISE survey gives, in addition to the magnitudes in the four 
bands, error estimates for these four magnitudes. To check the 
data as well as the error es timates, we compare d the w4 (22 yum) 
band data with data from iGautier et alJ (l2007h . who conducted 
a similar study using observations with the Spitzer space tele- 
scope. Some of our samples overlap (the corresponding stars are 
marked with C) in Table|2]). Spitzer, providing pointed observa- 
tions, performs significantly better than WISE in measuring in- 
frared magnitudes. For the 34 stars of which both Spitzer 24yum 
and WISE 22yum data ARE available, the Spitzer formal uncer- 
tainties (not i ncluding systematic calibration uncertainties, c.f. 
iGautier et all l2007h are consistently smaller, often by a factor 
exceeding 10. We can therefore use the Spitzer data as a bench- 
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Table 2. Observational data 
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/I Q/1 -uO OQ 


A A/l-uO 01 


/I '^/l -uO O^/ 
4. j4+U.Uj 


T 1 nm 
OJ 1 uuz 


1 4.oy±U.Uo 


Mj.j 




1 1 nn 
1 J. / / 


n _i_n no 

/ .44+U.UZ 


7 1 C-uO O^ 
/ . 1 o±U.Uj 


A CA-uO 00 


A 77 -uO OC/^ 
0. / / +U.UO-' 


Lrlo Z88 


1 4. / /+U.UD 


Mj.j 




1 1 oo 
1 j.VZ 


7 'TQ-i-n m 
/. / j+U.Uj 


7 '^O-i-O O^ 


7 OQ-uO 00 


A 7'^-l-0 (Ylf 

0. / J+U.O/-' 


T /I 1 O A 

OJ 4iZ A 


1 /I oA_i_n no 
1 4.oD±U.Uz 


A/T 1 A 
Ml.U 




a.l 1 


/I 'T'T.Ln no 
4. / /±U.Uz 


A A/1 _i_A AO 


A Af^^C\ Al 


/I A_i_n no f 
4. jD+O.Oz' 


OJ 412 D 


1 4.86+0.02 


M5.5 




\A.AA 


T O /I 1 A AO 

7.84+0.03 


T i 1 A AO 

7.61+0.02 


T 00 1 A AO 

/.zz+O.Oz 


T 1 1 A AA f 

7.12+0.09^ 


GJ 388 


1 A OC\ 1 A AT 

1 4.89+0.07 


M4.0 


4 


9.29 


A CA 1 A AO 

4.59+0.02 


y( y( O 1 A AA 

4.42+0.09 


^ , A AO 

4.27+0.02 


y( 1 C 1 A AO f 

4.15+0.03^ 


GJ 832 


1 4.95+0.02 


Ml. 5 


3 


8.66 


4.50+0.02 


4.29+0.08 


4.16+0.01 


4.05+0.02^ 


LP 944-020 


1 4.97+0.10 


M9.0 


2 


18.69 


9.55+0.02 


9.13+0.02 


8.27+0.02 


8.00+0.11^ 


GJ 166 C 


1 4.98+0.01 


M4.5 


3 


11.24 


5.96+0.03 


5.81+0.05 


5.48+0.02 


5.34+0.04^ 


GJ 873 


1 5.05+0.02 


M3.5 


1 


10.22 


5.30+0.02 


5.24+0.06 


4.89+0.02 


4.75+0.03 


GJ682 


1 5.05+0.05 


M3.5 


3 


10.95 


5.61+0.02 


5.35+0.07 


5.20+0.02 


5.13+0.04^ 


GJ 1116 A/B 


2 5.23+0.07 


M5.5 


1 


13.65 


6.89+0.02 


6.64+0.04 


6.28+0.02 


6.21+0.08 


G 099-049 


1 5.24+0.05 


M3.5 


1 


11.31 


6.04+0.02 


5.94+0.05 


5.61+0.02 


5.47+0.04 


LHS 1723 


1 5.32+0.04 


M4.0 


3 


12.22 


6.74+0.02 


6.53+0.05 


6.23+0.02 


6.13+0.05 


GJ445 


1 5.34+0.05 


M3.5 


1 


10.79 


5.95+0.03 


5.75+0.06 


5.51+0.02 


5.38+0.03 


GJ 526 


1 5.41+0.02 


M1.5 


1 


8.46 


4.42+0.02 


4.37+0.10 


4.19+0.01 


4.10+0.02 


GJ 169.1 A 


1 5.54+0.02 


M4.5 


5 


11.04 


5.72+0.02 


5.54+0.06 


5.32+0.01 


5.20+0.03 


GJ 251 


1 5.61+0.05 


M3.0 


1 


10.02 


5.28+0.02 


5.16+0.07 


4.92+0.02 


4.78+0.03 


2MA 1835+3259 


1 5.67+0.02 


M8.5 


6 


18.27 


9.17+0.02 


8.80+0.02 


8.16+0.02 


7.89+0.13 


GJ 205 


1 5.68+0.03 


M1.5 


1 


7.95 


4.04+0.26 


3.74+0.12 


3.70+0.02 


3.63+0.02 


LP 816-060 


1 5.71+0.11 


M4.0 


7 


11.50 


6.20+0.02 


6.02+0.05 


5.77+0.02 


5.63+0.04 


GJ 229 A 


1 5.75+0.03 


Ml.O 


1 


8.14 


4.17+0.23 


4.06+0.13 


3.97+0.02 


3.82+0.04 


GJ213 


1 5.83+0.03 


M4.0 


1 


11.57 


6.39+0.02 


6.23+0.05 


5.90+0.02 


5.68+0.04 


GJ693 


1 5.84+0.08 


M2.0 


3 


10.76 


6.02+0.02 


5.92+0.06 


5.60+0.02 


5.48+0.04 


GJ 752 A 


1 5.85+0.02 


M3.0 


1 


9.10 


4.67+0.02 


4.47+0.08 


4.38+0.01 


4.27+0.03 


GJ 752 B 


1 5.85+0.02 


M8.0 


1 


17.45 


8.77+0.02 


8.47+0.02 


8.08+0.02 


N/A 


GJ 570 B/C 


2 5.86+0.02 


Ml.O 


1 


8.09 


3.80+0.23 


4.06+0.06 


3.74+0.02 


3.62+0.03 


GJ754 


1 5.91+0.05 


M4.5 


3 


12.23 


6.85+0.03 


6.65+0.03 


6.36+0.02 


6.20+0.06 


GJ588 


1 5.93+0.04 


M2.5 


3 


9.31 


4.76+0.02 


4.70+0.07 


4.46+0.01 


4.36+0.03 


GJ 1005 A/B 


2 5.94+0.03 


M4.0 


1 


11.49 


6.39+0.02 


6.17+0.04 


5.91+0.02 


5.81+0.04 


GJ908 


1 5.95+0.04 


Ml.O 


1 


8.98 


5.04+0.02 


5.02+0.07 


4.76+0.02 


4.66+0.03 


GJ285 


1 5.98+0.07 


M4.0 


1 


11.19 


5.70+0.02 


5.51+0.06 


5.29+0.01 


5.13+0.03 
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Table 2. Observational data (contd.) 





stars in 


Distance 


Spectral 


SpTy 


V 


K,, 


^3.4 


Fl2 


F22 


Name 


beam 


(pc)° 


Tvpe* 


ref. 


(mag)'^ 


(mag)'' 


(mag)"^ 


(mag)'' 


(mag)'' 


GJ 268 A/B 


2 


6.12±0.07 


M4.5 


1 


11.48 


5.85+0.02 


5.69+0.06 


5.38+0.02 


5.20+0.03 


GJ784 


1 


6.20±0.04 


K7.0 


3 


7.95 


4.28+0.02 


4.15+0.10 


4.05+0.02 


3.98+0.02 


GJ555 


1 


6.22±0.08 


M3.5 


1 


11.32 


5.94+0.03 


5.79+0.06 


5.54+0.02 


5.41+0.03 


GJ 896 A/B 


2 


6.25+0.06 


M3.5 / M4.5 


1 


10.10 


4.94+0.03 


4.76+0.08 


4.64+0.02 


4.59+0.03 


GJ581 


1 


6.34±0.07 


M2.5 


1 


10.57 


5.84+0.02 


5.69+0.06 


5.48+0.01 


5.33+0.03 


LHS 2090 


1 


6.37±0.11 


M6.0 


2 


16.10 


8.44+0.02 


8.23+0.03 


7.80+0.02 


7.96+0.22 


LHS 337 


1 


6.38+0.08 


M4.5 


3 


12.75 


7.39+0.02 


7.24+0.03 


6.86+0.02 


6.72+0.05 


GJ 661 A/B 


2 


6.40+0.05 


M3.0 


1 


9.37 


4.83+0.02 


4.71+0.08 


4.46+0.02 


4.34+0.02 


G 180-060 


1 


6.41+0.16 


M5.0 


3 


14.76 


8.51+0.02 


8.29+0.02 


7.95+0.02 


7.90+0.13 


GJ 644 A/B/D 


3 


6.45+0.02 


M2.5 


1 


9.03 


4.40+0.04 


4.22+0.09 


4.04+0.02 


3.90+0.03 


GJ 644 C 


1 


6.45+0.02 


M7.0 


1 


16.78 


8.82+0.02 


8.59+0.02 


8.13+0.02 


7.86+0.18^ 


GJ643 


1 


6.45+0.02 


M3.5 


1 


11.80 


6.72+0.02 


6.57+0.03 


6.29+0.02 


6.12+0.05 


GJ625 


1 


6.53+0.04 


M1.5 


1 


10.10 


5.83+0.02 


5.68+0.06 


5.46+0.02 


5.33+0.02 


GJ 1128 


1 


6.53+0.10 


M4.5 


1 


12.74 


7.04+0.02 


6.80+0.03 


6.51+0.01 


6.40+0.03 


ni 1 1 S6 

VJ J 1 i JU 


I 




M5.0 


]^ 


lJ.O\J 


7 S7+n fl^ 

/ .J> / ZIZKJ.KJJ 






\j.y 1 mu.uo 


LHS 3003 


1 


6.56+0.09 


M7.0 


1 


17.05 


8.93+0.03 


8.69+0.02 


8.27+0.02 


8.12+0.27^ 


GJ408 


1 


6.70+0.07 


M2.5 


3 


10.02 


5.50+10.00 


5.38+0.06 


5.19+0.01 


5.09+0.03 


GJ 829 A/B 


2 


6.71+0.08 


M3.5 


1 


10.30 


5.45+0.02 


5.30+0.07 


5.06+0.01 


4.90+0.03 


G 041-014 A/B/C 


3 


6.77+0.09 


M3.5 


3 


11.46 


5.69+0.02 


5.49+0.06 


5.25+0.02 


5.08+0.03 


GJ402 


1 


6.80+0.14 


M4.0 


1 


11.65 


6.37+0.02 


6.26+0.05 


5.98+0.02 


5.86+0.04 


GJ880 


1 


6.83+0.04 


M1.5 


1 


8.65 


4.52+0.02 


4.43+0.08 


4.33+0.01 


4.22+0.02 


GJ299 


1 


6.84+0.14 


M4.0 


1 


12.82 


7.66+0.03 


7.44+0.03 


7.12+0.02 


7.01+0.09 



Notes. 

Distances calculated from RECONS parallaxes iHenrv et al.Ll200d) . 

Spectral types from the literature. When only one spectral type is given for a binary / multiple system, this represents the spectr al type of the 
whole s ystem and the indiv idu al spectral types are un known. The sources of the sp ectral type are giv en in the next column a nd are: DiHenrv et all 
I200I 2) lHenrv et allllOOl 3) iHawlev et .511199^, 4) lKeenan & McNeill [1981 5) lGraY et al.Ll2003l. 6) lReid etalll200l 7) lGrav et alll200l 
'''' V magnitudes talcen from the RECONS database; where more than one star is in the WISE / 2MASS beam, the fluxes are added. 

A", magnitudes and errors are taken from the 2MASS database dSkrutskie et aU[2006h ; where more than one star is in the WISE beam, the fluxes 

of all stars in the beam are added. 

''"' WISE bands wl centered at 3.4yum, w3 centered a t 12/jm, and w4 cent ered at 22^m jWright et al.Ll2010l) . 

Stars with 24 yum Spitzer flux measurements from lGautier et alj (l2007l) . 



mark for the WISE data, as we did in Figure |2] There is a slight 
offset between the Spitzer and WISE magnitudes that we deter- 
mined to be 0.083 magnitudes using a weighted least-squares 
fit. For the high-luminosity end, the fit is very good and suggests 
that the WISE errors might actually be overestimated (or there is 
a fraction in these errors that represents the global calibration un- 
certainty, which has no influence on this comparison other than 
a global offset, which we corrected for). At the low-luminosity 
end, there are some outliers. At least one of these outliers can be 
explained because visual inspection shows that the WISE source 
is not exactly centered on the actual source. We conclude from 
this that at least in the longest waveband, the data quality is 
mostly good, but some outliers might be expected. No such state- 
ment can be made about the other wave bands. However, as our 
analysis will show, the data quality is indeed excellent and the 
WISE errors seem to be generally overestimated. 



3. Analysis 

For our 85 systems, we obtained magnitudes in seven different 
bands (V, the 2MASS bands J, H and K, as well as the WISE 
bands wl, w3 and w4). In addition, we also have parallax and 
thus distance information, giving us access to the absolute mag- 
nitudes. The difficulty we face with M-dwarfs is that, unlike 
solar-type or A-stars, they are not well into their Rayleigh-Jeans 
tail at 3.4 jum. This means that not only the J-H, J-Ks and H-Ks 



colors depend on spectral type, but that the same applies for the 
colors in the WISE bands. 



3.1. A proxy for stellar temperature 

Because we expect the dependence of the WISE colors for stars 
without excess to be mainly with the effective stellar temperature 
Tgff, we need a good proxy for this. Casagrande et al. (2008!) 
have shown that none of the 2MASS colors (J-H, J-Kg, H-K,) 
are good proxies for effective stellar temperature, but show sig- 
nificant scatter and are not necessarily a monotonic function of 
spectral type. Good proxies are the V-Ks color as well as abso- 
lute magnitudes in a variety of bands. 

The spectra of M-dwarfs are dominated in the visible TiO 
bands. The V-band flux for a given 7^// may therefore sig- 
nificantly depend on abundances and atmospheric dynamics. 
Furthermore, man y M-dwarfs are known to be variable in the 
visible band (e.g.. Ha wlevl Il993f ). which in turn can influence 
the V-Ks color, depending on at which point in the variability 
cycle the V-band data was taken. To assess this, we plot the ab- 
solute 3.4 fim magnitude obtained from WISE against the V-Ks 
color in Figure[3] If both are indeed good proxies for stellar tem- 
perature, we expect to be able to fit a (not necessarily straight) 
line through the data points. This is clearly not the case. Plots 
against the K, and V absolute magnitude show that the main 
contribution to the scatter comes from the V magnitude. 
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Fig. 3. Comparison of WISE absolute magnitudes in wl (3.4 /zm) band with V-Ks color (left panel). As can be seen, the two proxies 
for stellar temperature do not match well. This is not because of the Ks magnitude (as can be seen in the middle panel), but because 
of the magnitude in the visual band (right panel). An empirical determination of M-dwarf mid-infrared colors as a function of the 
V-Ks color index can be found in the appendix. 



We conclude that because of strong molecular absorptions 
and variability in the V band, the V-Ks color is not an ideal 
proxy for the stellar temperature of late-type dwarfs. We chose 
to use the 3.4 fj.m absolute magnitude instead, because it is a 
good proxy for stellar temperature and furthermore can give us 
the best data consistency because it originates from the same in- 
strument. The reduced goodness-of-fit measurements of our 
analysis support this choice because the fits to the data are better. 

There is one caveat to using absolute magnitudes as a proxy, 
though. Because pre-main-sequence stars are too bright for their 
spectral type due to their larger diameters, the predictions for 
their spectral types / surface temperatures and consequently their 
colors in the mid-infrared will be off (c.f. Section [TSl l. For this 
paper, this is not a problem because there are no known pre- 
main-sequence stars in our sample. For young stars, however, it 
is better to use the V-Ks color as a proxy, which is indeed also 
a very decent proxy for stellar temperature even for these late- 
type dwarfs. We show the fit against the V-Ks color index for 
comparison and reference in the appendix. 

3.2. Variability in fhe WISE bands 

As we have discussed above, many M-dwarfs are known to be 
variable in the visible. Several of our stars are flare stars, produc- 
ing short-term variability in the ultraviolet. Out of t he 34 stars in 
our sam ple that were also observed with Spitzer bv .Gautier et alJ 
(l2007h . the SIMBAD database marks more than half as flare 
and/or variable stars. While we cannot comment on the accuracy 
of these classifications, we see that flare and variable stars are 
part of this sample. On the other hand, we see that the WISE and 
Spitzer magnitudes match well (cf. Figure |2]i at the bright end , 
while outliers are present at the faint end. Toffl emire et al.l(l2012l) 
observed three M-dwarfs for flares and short-term variability and 
found that while the flares are intense in the ultraviolet, the vari- 
ability they produce in the near-infrared 2MASS regime is be- 
low 7-11 mmag depending on the band. From this we conclude 
that the outliers in the longer WISE bands are mostly not due 
to variability, but have other reasons (for example, bad center- 
ing, background contaminations or contaminations from nearby 
sources). 

We can get rid of most of the variability in the mid-infrared 
by using the WISE colors. All four WISE bands are taken at the 
same time, which means that any variability that influences all 
mid-infrared bands equally will not influence the WISE colors at 



all. The WISE catalog provides a variability flag for each of the 
four bands. Out of the 85 sources in our sample, only a small 
minority of seven are flagged as variable in at least one of the 
bands. We conclude that variability is likely to not be an issue 
for excess detection in our analysis for most of our sample. 

3.3. Multiplicity 

As we have noted before, WISE is not able to resolve all known 
multiple systems within our sample. Because we do not want 
to exclude these systems from our study, we have to investigate 
what happens to the WISE magnitudes and colors in multiple 
systems. 

Equal-mass binaries will simply have double the luminos- 
ity of the single star, the colors should not be influenced. They 
move to the left by 0.75 magnitudes in FigureH] Binaries of dif- 
ferent mass will move less in absolute magnitude (compared to 
their primary alone), but might move slightly upwards if the sec- 
ondary contributes significantly and has a significantly different 
color This effect gets smaller for a larger difference in spectral 
class and thus a larger color difference because the secondary 
will then only have a fraction of the brightness of the primary, 
and we expect these binaries to lie relatively close to the posi- 
tion their primary would occupy in the diagram. The same is in 
principle true for triple systems, even though somewhat larger 
deviations are to be expected. 

Because of the above discussion, we do not use the multiple 
systems within our sample to obtain the polynomial fit. We do 
compare them to the fit, though - bearing in mind that we will 
have to treat them specially if any of them turns out to show 
signs for excess. Binary and triple systems are marked in Figure 
|4] As we expect, the multiple systems lie preferentially above 
the fitted polynomial line. 

3.4. Fitting and intrinsic spread of mid-infrared colors 

In order to be able to detect excesses, we need to know the 
expected mid-infrared color for every source. Instead of using 
model stellar atmospheres, we choose a purely data-driven ap- 
proach and use the data at hand (which might include excess 
sources) to fit for the expected mid-infrared colors. This is only 
possible if the number of data points significantly exceeds the 
number of parameters to fit. It also assumes that the number 
of excess sources is small compared to the general sample size. 
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Fig. 4. WISE colors plotted against 3.4 fim absolute magnitude. Single M-dwarfs are blue, binary systems are marked in red, triple 
systems in green. The black open triangles show the K dwarfs included to better constrain the high-mass end of the fit. The black 
dashed line shows the third-order polynomial fit using all single dwarfs, while the black solid line shows the fit after iteratively 
removing outliers and is used for our analysis. The inset shows the histogram for the cr' error value adopted for our analysis (see 
definition in the text). The solid bars represent the stars that were used in the fit, the blank stacked bars are for the stars that were 
omitted from the fit (multiples and outliers). In all three cases, the histograms are consistent with a standard-normal distribution 
(plotted here as a red line). For reference and as test of our technique, the black filled downward triangle to the left represents the 
well-known excess source AU Mic, which seems to have significant excess in all three colors. For a discussion of AU Mic, see 
SectionllSl 
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Table 3. Statistical tests for cr' 



Color index 


KS test 


SW test 




(p value) 


(p value) 


[3.4]-[12] 


0.84 


0.82 


[3.4]-[22] 


0.74 


0.71 


[12]-[22] 


0.47 


0.41 



Notes. Probability of the data (cr') being drawn from a Gaussian dis- 
tribution as calculated using the KS (Kolmogorov-Smimov) and SW 
(Shapiro- Wilk) tests. All data are consistent with being standard-normal 
distributed. 

which is a reasonable assumption (c.f. section|4]i. We use a poly- 
nomial of 4th order to describe the mid-infrared colors and ob- 
tain the fit through an iterative fitting routine that removes out- 
liers in each step (details on this fitting routine and its parameters 
can be found in the appendix). Outliers are defined as sources 
that diverge from the fit by more than 3cr'. cr' is defined from the 
data and the fit as 



where cr, is the WISE error estimate for source i, e, is the devi- 
ation from the fit and is the overall goodness-of-fit measure 
of the fit. We note at this point that the iterative fitting proce- 
dure does not change the number of excess detections (i.e. there 
is no excess source that would not have been detected as excess 
so urce had we no t used an iterative fit). 

iGautier et alJ (2007) stated that the limits for detecting mid- 
infrared excesses stem from the scatter in photospheric flux 
predictions from shorter wavelengths and that the mid-infrared 
properties of M-dwarfs are not well known. While we did not 
predict the mid-infrared fluxes, but instead directly fit the col- 
ors, we are also affected by a possible intrinsic spread in the 
mid-infrared colors, the amount of which is unknown. However, 
we know that if the intrinsic spread was negligible, our model 
was a reasonably good fit to the data (which visually seems the 
case) and the WISE error estimates were exactly right, then the 
deviations from the fit divided by their respective WISE error 
estimates (|f-) should follow a standard-normal distribution and 

the measure should be ~ 1 . In fact, the measures we ob- 
tained were significantly smaller than one, which means that the 
intrinsic spread must be significantly smaller than the WISE er- 
ror estimates, and furthermore, the WISE errors are presumably 
overestimated. 

Because separating the errors from the (overestimated) WISE 
measurements and the intrinsic spread would only be possible 
with a large set of (not necessarily justified) assumptions, we 
used the above-defined cr' value instead. If the amount of overes- 
timation is about the same for all WISE error estimates and if the 
intrinsic spread is in fact much smaller than these estimates, cr' 
follows a standard-normal distribution. We present histograms of 
cr' in Figure|4]as insets. We also used the Kolmogorov-Smirnov 
test and the Shapiro- Wilk test to check whether the distribu- 
tions are compatible with a Gaussian distribution. The results 
are shown in Table [3] For all three examined colors, there are no 
hints that the distribution is not Gaussian. 

It is worth noting that we fit the three colors individually and 
not simultaneously, which means that they are not self-consistent 
and the net color ([3]-[12] + [12]-[22] - [3]-[22]) is not zero. The 
residuals, however, are so low (0.005 magnitudes on average) 
that this has no influence on our excess detection technique. 



3.5. Excess detection 

Because of the Gaussian nature of cr', it can directly be used 
as a measure to detect excess. We considered a source to show 
statistically significant excess if cr'. > 3. No source in any of the 
three colors shows a statistically significant deficit, defined as 
cr'j < -3. This gave us confidence that significant excess would 
not be just a statistical outlier 

3.6. Excess sources 

In our sample, we found two sources with cr' values > 3 in [3.4]- 
[12] color, three in [3. 4] -[22] color and one in [12] -[22] color. 
We list these sources in Table |4] The six detections are for four 
individual stars. We discuss each star in detail. 

3.6.1. GJ 570 B/G 

The binary GJ 570 B/C shows excess at 12 and 22//m and is 
flagged as variable in all three considered bands. Visual inspec- 
tion of the WISE images shows a second nearby source. The GJ 
570 system actually consists of four individual stars: The K4 pri- 
mary GJ 570 A, the two close M-dwarfs GJ 570 B and C, and 
the T7 brown dwarf GJ 570 D. The brown dwarf is at large sep- 
aration (~ 250"), the source next to the B/C binary is the A pri- 
mary at a separation of about 20". Both sources are bright (3.8 
magnitudes in Kg), and since they are close, contamination could 
be a problem. Indeed, a close examination of the J, H and Ks 
band in combination with the WISE bands yields that the GJ 570 
B/C is measured to be more than 0.25 magnitudes fainter in the 
first WISE band than in Ks and more than 0.3 magnitudes fainter 
in the first WISE band compared to the third. We conclude that 
the most likely explanation is that the measurement in the first 
WISE band is incorrect. This explains both the apparent excess 
in the third and in the fourth band, since the colors are calculated 
against the first WISE band. It would not explain any excess in 
[12] -[22], and indeed no excess is seen at this color (the excess 
value of 1.7cr' is well below our threshold of 3cr' and might be 
artificially enhanced because GJ 570 B/C is a binary, c.f. section 
13.31 ). We therefore conclude that there is no significant sign for 
excess at either 12 or 22jum. 

3.6.2. SCR 1845-6357 A 

This single star shows a clear 4.8cr' excess at 12//m and no ex- 
cess at 22jL/m, although the associated cr' value is 1.9. The de- 
rived cr' of 0.87 in [12]-[22] is not significant. It is not flagged 
as variable in any of the WISE bands. The excess in magnitudes 
over the fitted line is 0.132, 0.144 and 0.064 magnitudes, corre- 
sponding to 12.9%, 14.2% and 6.1% excess, respectively, with 
the caution that the latter two are not significant. 

The excess can be explained by the fact th at SCR 1845-6357 
has a T6 secondary at a separation of ~ 1 " dBiller et al.L 120061 : 
iKasper et al.L |2007|) . While T dwarfs usually do not contami- 
nate the flux of a system sufficiently to produce a detectable ex- 
cess in solar-type stars because of the large difference in emitted 
flux, S CR 1845 -6357 A is a very late M-dwarf of spectral class 
M8.5 (iHenrvetalll200 6) and thus itself very faint . The T dwarf 
is fai nter by only ~4 magnitudes in the H band dKasper et all 
i2007h . While the flux in the 3.4//m band of WISE is only on flie 
level of ~ 1 - 2% of the M-dwarf primary, it can explain the 
increased flux in the 12jL/m and 22yum bands. We conclude that 
the excess is real in both the 12 and 22jt/m bands, but does likely 
not originate from a debris disk, but from the T dwarf secondary. 
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Table 4. Sources showing excess in our analysis - discussion in text 



Name 


[3.4]-[12] 


[3.4]-[22] 


[12]-[22] 


Comment 


GJ 570 B/C 


5.27 


6.05 


1.73 


Contamination from nearby K4 primary 


SCR 1845-6357 A 


4.84 


1.92 


0.87 


T6 brown dwarf companion in beam 


LHS 288 


1.34 


3.78 


3.49 


22yum measurement strongly inconsistent with Spitzer measurement 


GJ 860 A/B 


2.41 


3.21 


2.35 


Apparent excess caused by binarity 



Notes. The four sources identified as excess sources from our tr' > 3 test with their corresponding cr' values for the three colors we constructed. 
Statistically significant excesses are marked in bold. The comment gives possible reasons why the detected excess could be spurious or not originate 
from a debris disk. The sources are discussed individually in the text. 



This shows the validity of our method, which is indeed able to 
pick up excesses in the 12//m band of as low as 13% for this very 
late dwarf. 



3.6.3. LHS 288 



LHS 288 is a single star showing excess in [3.4]-[22] (3.8 cr') 
and [12]-[22] (3.5 cr'). There is no clear indication that the 22jt/m 
measurement of 6.75 magnitudes would be wrong, even though 
it can be seen in the WISE image that there is a diffuse nebulos- 
ity in the background. The star is furthermore flagged as variable 
in the 3.4 and 22 fim WISE bands. It is worth noting, however, 
that th is source was also measured with Spitzer bv lGautier et alJ 
(I2007h . The two measurements are strongly inconsistent, with 
the WISE measurement being 9.5 cr above the Spitzer measure- 
ment at 24jL/m. We conclude that the most likely explanation for 
LHS 288 is a poor measurement in the last WISE band and not a 
real excess. 



3.6.4. GJ 860 A/B 



The close M3/M4 binary GJ 860 shows a derived excess at the 
longest WISE band of 3.2cr'. This corresponds to 0.254 magni- 
tudes or 26.3% excess at this wavelength. At 12pm, the excess is 
at a 2.6cr' level (corresponding to 0.184 magnitudes / 18.4%), 
the signal in [12]-[22] color is 2.3cr' corresponding to 0.067 
magnitudes or 6.4%. The WISE catalog flags this source as vari- 
able, but only in the third (12yum) band. The system does lie in a 
crowded region with possible background cirrus, but it is bright 
itself (magnitude 4.7 at 3Apm, magnitude 4.1 at 22pm). Other 
bright sources are nearby, especially one source that is strong 
at 22pm but weak at the other wavelengths (in comparison to 
GJ 860). This source lies approximately 90" to the southwest. 
It is unlikely that it strongly affects the WISE measurement at 
this wavelength. Because this is a binary star, special precau- 
tions are necessary. The increased luminosity of the total system 
will move the star to the left in our diagram, and because the 
later-type member of the system will show a more red [3.4]- 
[22] color, it will shift slightly upward compared to a single M3 
dwarf. Assuming that the M3 dwarf is stronger by a factor of 2 
in total luminosity, both effects can be calculated and the posi- 
tion of the M3 dwarf alone in the diagram can be derived. Doing 
so, we derive new values for cr' of 2.3, 1.7 and 1.9, respectively. 
None of these are significant, and we therefore did not include 
this source in our excess sources. An excess might be present at 
levels below our detection threshold. 



3. 7. Excess detection limits 

Because we handled each star individually rather than setting a 
general limit on F/r/F*, with Fjr being the detected infrared 
flux and the expected stellar flux from our model fit, we did 
not obtain a global value of an excess detection limit for all our 
stars. Instead, we were able to detect (3cr' limit) infrared flux 
excesses down to 10% in some stars while in other stars, even a 
30% excess would not be detected due to a larger WISE measure- 
ment error. Because of this, we calculated the excess detection 
limit for each star individually for each color assessed from its 
WISE error and the value of the corresponding color fit. The 
typical (median) limits we derived are 14.9% for the [3.4]-[12] 
excess, 17.9% for the [3.4]-[22] excess and 7.0% for the [12]- 
[22] excess. The mean values are higher due to a few sources 
for which we derived quite poor detection limits and measure 
15.4%, 20.2% and 10.9%, respectively. 

These limits can be converted into limits on the fractional 
dust luminosity L/s/L* if we assume a single dust ring of a given 
temperature (see Fig.|5]l. As can be seen, we put no strong con- 
straints on the fractional dust luminosity for cold dust of tem- 
peratures lower than ~50K. Large amounts of cold dust could be 
present around the M-dwarfs we observed to have no significant 
amount of excess at 12 and 22pm. At our most sensitive point, 
which is in the interesting (because of its link to terrestrial planet 
formation zones) 100-300K range, we derived 3cr upper limits 
for L/r/L* of nearly 10"'* for Barnard's star (M3.5, 1.8pc) and 
nearly 6 * 10"^ for GJ 526 (Ml, 5.4pc). These are indeed fairly 
typical values for our survey, although they depend slightly on 
the spectral type - for the earliest stars in our sample, we can go 
deeper by a factor of ~2, while for the latest stars we are less 
deep, again by a factor of ~2 (both compared to Barnard's star). 
This is because it becomes increasingly easy to detect dust when 
the differ ence in temperature betwee n the dust and the host star 
is larger (iBackman & P ares^. '1 9931). Also be cause of this, fol- 
lowing the discussion in Gautie r et al.l (l2007l) . we realized that 
while the limits we obtained for the fractional dust luminosity 
are less stringent than have been derived for earlier-type stars, 
this indeed corresponds to more stringent constraints on the dust 
mass at these temperatures (see discussion in section l43T l. 

In most of the cases, the [12] -[22] color is in the most sensi- 
tive to dust of these temperatures because such dust produces no 
significant excess at I2pm and because our excess limits are the 
most stringent ones in this color. 

3.8. AU Mic 

Although not in our sample due to its distance, it is interesting 
to investigate in passing the most well-known M-dwarf with a 
debris disk, AU Mic, an ~12 Myr old Mle star in the /? Pictoris 
moving group dZuckermani 1200 ih . While excess emission has 
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- 3a limit from 12 \im excess 
3a limit from 22 fim excess 

- 3a limit from 22 to 1 2 nm color 




10" 



Dust temperature (K) 




Dust temperature (K) 

Fig. 5. Inferred dust limits for GJ 699 (Barnard's star, top) 
and GJ 526 (bottom). The different lines show the limits from 
the different WISE colors. The gray shaded region is excluded 
by the combination of all three individual limits. For refer- 
ence, the black dashed line shows the expected performance of 
SPHERE/ZIMPOL for a thin dust ring (see section|5]l. 



been found at wavelengths of IQ/im and beyond, no excess has 
been rep orted at mid- infrared wavelengths of 24/vm or shorter 
jLiu et al.. 2004; Plavchan et al., 2009). We plot it together with 
the stars from our sample in Figure |4] where it can be seen to 
have significant excess in all colors (3. Oct', 4.9cr' and 5.8cr', re- 
spectively). However, it lies in the K dwarf regime w.rt. its ab- 
solute 3.4yum magnitude. This is because it has a larger surface 
area than the average Ml dwarf due to its young age. As we have 
seen, the M-dwarf mid-infrared colors depend on their spectral 
type / surface temperature, and the absolute 3.4jum magnitude is 
a poor proxy for either in case of pre-main-sequence stars that 
are too bright for their spectral type. 

It is accordingly better to use the V-Ks color as the proxy 
for stellar temperature for such young stars. Doing so, perform- 
ing the fits and comparing the colors of AU Mic to them, we 
obtained cr' values of 0.3, 1.2 and 3.3, respectively, for the [3.4]- 
[12], [3.4]-[22] and [12]-[22] colors. This is significant in the last 
color, pointing to an excess of 6.7+2.0% over the photosphere at 



22/j.m. This can be compared to lPlavchan et alJ (l2009l) . who did 
not find a significant excess (3.5% above the photosphere, 1.65cr 
using Spitzer MIPS) and Simon et al. (2012), who did not detect 
an excess at 22/im either using the WISE preliminary release cat- 
alog and a fixed threshold of [12]-[22]>1.0, stating that the AU 
Mic debris disk is too cold for detection with WISE. 

This means that no mid-infrared excess is known for AU 
Mic, but both these results are consistent with its disk hav- 
ing an excess at 22/24jum, albeit a small one. The Spitzer 
MIPS measurement of 155.2+3.2 mJy can be converted to 
4.160+0.028 magnitudes using the zero points from the Spitzer 
MIPS Instrument Handbook. This is consistent with the WISE 
measurement of 4.137+0.025 magnitudes, both with and with- 
out taking into account the general offset between the Spitzer 
and WISE measurements discussed in Section l2.2.1l This means 
that we only detect the small excess because our method is more 
sensitive. Indeed, comparing the two bands of WISE directly us- 
ing a calibration from other stars might currently be the most 
sensitive way to detect very small excesses. 

Additional hints of an excess at shorter wavelengths than 
previously thought for AU Mic come from an unpublished 
Spitzer IRS high-resolution spectrum of AU Mic (5-35;um) taken 
in 2005. That spectrum has unfortunately never been published, 
and the data reduction and especially background subtraction is 
not trivial, but a quick glance at that spectrum does hint at an 
excess starting at ~18/im. Unfortunately, the IRS spectrum and 
the MIPS measurement are not compatible, and resolving the 
apparent contradiction is difficult. 

There are more caveats to this. Looking at the images in the 
WISE image archive, the centering of the source seems to be in- 
accurate. Because of this, are cautious about these measurements 
of AU Mic. AT Mic, another M4.5 member of the yS Pictoris 
moving group, detected with a n excess of ~15% at lAjim using 
Spitzer ( Plavchan et al.L 120091) . is not detected as having excess 
using our technique. It is also not perfectly centered in WISE. 
It has an absolute 3.4/im magnitude similar to that of AU Mic, 
pointing out again that absolute magnitudes are a poor proxy for 
stellar temperatures / spectral types of pre-main sequence M- 
dwarfs. 

4. Comparison to earlier studies 

Several studies have been performed on infrared excesses 
for main-sequence stars. We concentrate here on studies of 
25//m and shorter-wavelength excess. The first major sur- 
vey of this kind was conducted with IRAS. However, IRAS 
was on the whole not s u fficiently sensit ive for M-dwarfs 
(iBackman & Parescel Il993h . ISu et al.l (12006 ') searched -160 A- 
stars for excess w ith Spitzer, yielding an excess rate of 32^^%. 
iMoor et aP (|201 ih concentrated on F-type stars and observed 82 
stars with Spitzer, of which 15 show excess at 24yum (18^^%), 
but in this case, the sample selection was biased toward stars 
that had shown hints for excess in other surveys. A more re- 
cent survey combining Spitzer and WISE data for a total sam- 
ple of 263 F-stars detected excess in 19 sources , correspond- 
ing to 7.2+[-^% (Mizusawa et al., 2012, submittedl. lTrilling et all 
(2008) observed nearly 200 sun-like (FGK) stars and found that 
the incidence of debris disks detectable at 24/zm is 4.2^j-^%. 
These authors also specifically targeted binary and multiple 
early-type (A3-F8) systems and found that the incidence of de- 
bris disks is comparable to or slightly higher than for single stars 
(Trilling et al., 2007). Carpenter et al. (2009) also targeted solar- 
type stars, concentrating on the evolution of detectable excesses 
with age. They found that the fraction of systems with detectable 
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excesses (> 10.2% above the photosphere) at 24//m drops from 
15+^% for ages < 300 Myr to only 2.7+|-^% for older ages. 
iBrv den et al. (2009) used Spitzer to determine whether the inci- 
dence of detectable excesses is different for stars known to har- 
bor planets from radial-velocity searches. They found that any 
observed difference to stars not known to harbor planets is not 
statistically significant, but add to the general picture by con- 
tributing observations of 104 stars of spectral types F through 
M of which four (2.7;^[ 2%) show detectable excesses at 24jum. 
None of their excess sources is an M-dwarf. 

All these s urveys targeted stars earlier than spectral type M. 
iPlavchan et al.l ([2005) pointed out that there is a lack of debris 
disk detections for mature M-dwarfs and discussed possible rea- 
sons. iGautier et al.l (l2007h performed a Spitzer search around 
nearby 62 M-dwarfs (some of their targets coincide with ours) 
and found no exc esses at 24/^m. 

The study by Fuiiw ara et al.l 12009) was performed using the 
data from the all-sky AKARI satellite. Their detection thresh- 
old was significantly higher at 50% above the photosphere (at 
18//m), where they found significant excess around 14 main- 
sequence stars and derived an excess fraction of 1.5%. None of 
their excess stars is a K or M-dwarf. 



4.1. Dust at similar temperature 

As described above, we detected no significant excess that holds 
up to critical assessment in any of our stars. Because we studied 
85 targets in the 3.4 and 12yum bands and 84 targets in the 22/im 
band, we derive an excess fraction (at the described levels) of 
0;^()Q% amongst M-dwarfs. We did not consider the physically 
real excess around SCR 1845-6357 A because it stems not from 
a debris disk, but from a T dwarf binary. T dwarf binaries will 
not produce detectable amounts of excess around AFGK-stars, 
given current photometric precisions. 

The discussion is somewhat complicated because for some 
of the targets, we have more than one target in the beam. 
Considering only targets with single stars in the beam (70 at 
3.4 and 12jum, 69 at 22pm) yields an excess fraction of 0^^^%, 
with some of them being part of wide-se paration (dis tinguish- 
able by WISE) binaries. However, Trillin g et alj (l2007h showed 
that the excess fraction in binary systems is comparable to or 
slightly higher than that for single stars for a sample of A and 
F main-sequence stars. Furthermore, not all studies mentioned 
above considered only single stars. Because of this, we decided 
to directly compare our derived excess fraction of O^g g% to the 
studies of AFGK and M-dwarfs performed in the literature. 

Our derived detection rate for M-dwarfs is lower than for 
any other spectral class. This effect cannot simply be explained 
by the fact that our survey is less deep. While we reach only a 
median excess detection limit of 17.9% in the [3.4]-[22] color (to 
be compared with usually ~ 10-15% for the Spitzer surveys), 
our median excess detection limit in the [12]-[22] color is 7.0%. 
In the most sensitive temperature range, the [12] -[22] color usu- 
ally gives the strictest limits (c.f. Figure|5]l. We conclude that we 
should have been able to detect a comparable percentage of ex- 
cess systems as would have been possible with a Spitzer survey 
of the same targets. 

Comparing our findings to the A-star survey of ISu et alJ 

(l2006l) . we can conclude that (randomly selected) M-dwarfs 
have a lower fraction of excess systems at a ~ 6cr level. 
Comparing our results to the F-star sample of Mizusawa et al. 
(2012, submitted), the difference of the combined M-dwarf sur- 
veys is significant on a ~ 3.6cr level. However, as shown by 
ICarpenter et al.l (l2009h . the fraction of excess systems is strongly 



dependent on age. A- and F-stars are preferentially younger (be- 
cause they live for a shorter time). Our sample of nearby M- 
dwarfs is likely to be sampled from the average population of 
M-dwarfs in our galaxy, we assume an average age of several 
Gyr (the star forma t ion ra te of the Milky Way peaks at ages 
-2-6 Gv r. seelWvsel l2009l). The stellar age of the sample from 
ISu et all (12006')" is on average only a few hundred Myr, with no 
star older than 1 Gyr. Mizusawa et al. (2012) provided no ex- 
act ages but estimated that their s ample is probably on average 
~ 300 Myr. 'Carp enter et all (12009 ) found that the excess frac- 
tion for FGK-stars drops to 3.6% for ages 300-1000 Myr and 
1.8% for ages >1000 Myr They also found that the upper enve- 
lope for 24yum excess at ages > 500 Myr is ~ 10%. If the same 
was true for our M-dwarf sample, of which the vast majority is 
expected to be older than 300 Myr and a majority to be older 
than 1000 Myr, this means three things: 1) The declining excess 
fraction with age can explain the scarcity of M-dwarf excess de- 
tections at 24pm and shorter wavelengths, 2) the excess fraction 
of M-dwarfs is not, on a 3cr level, statistically different from the 
excess fraction of FGK-stars, and 3) we might be missing some 
debris disks very close to our detection threshold. 

Our results agree with lGautier et al] (l2()()7i) . who did not find 
any excess around M-dwarfs. When combining the two samples 
for a total number of 1 12 M-dwarfs probed for excess at either 
22 or 24pm and making the conservative assumption that 25% of 
these M-dwarfs are in the 300-1000 Myr age range and the rest 
is > 1000 Myr old, the excess fraction of the M-dwarf sample 
is different from the FGK sample on a ~ 1.1 cr level, there is no 
clear evidence for a fundamentally different excess behavior of 
M-dwarfs at these wavelengths (corresponding to warm dust), 
but because of the general scarcity of M-dwarf excesses, more 
data are needed to clarify this question. 

4.2. Dust at similar radii 

In the last section, we have compared our 22yum (non)detections 
for M-dwarfs with studies of earlier-type stars at 24pm. This 
means, however, that we compared dust at very different radii. 
The orbital radius of dust of a given temperature, assuming 
blackbody grains, scales as Rdust ~ L^^^- This means that while 
we probed dust at scales <K 1 AU for our M-dwarf sample, the 
same wavelength corresponds to dust at ~ 1 AU for solar-type 
stars and > 1 AU for A-stars. 

The discussion is additionally complicated because the dust 
behaves distinctly different from blackbody grains. For M- 
dwarfs, blowout due to radiation pressure is not relevant because 
of their low luminosities. This means that the dominant grain 
removal proc esses will b e Poynting-Robertson (P-R) and stellar 
wind drag (Plavchan et al.. 2005: Saiia et al.. 2003). Depending 
on the stellar wind from the M-dwarfs, the stellar wind drag 
can dominate. Since the stellar winds of M-dwarfs are n ot well 
known, their effects are difficult to quantify (Plavchan et all 
2005) and the expected grain size distribution around M-dwarfs 
is unknown. On the other hand, grains smaller than the peak 
wavelength of their thermal emission will not behave like black- 
body grains, but instead have higher effective temperatures than 
blackbody grains at the same location around the star. This 
means that we are in fact probing orbital radii larger than ex- 
pected for blackbody grains with our observations (but likely 
still significantly smaller than the orbital radii we probe with ob- 
servations at the same wavelengths around earlier-type stars). 

Excess emission around stars of any type becomes rarer 
for shorter wavelengths. Excess at wavelengths shorter than 
22 pm around main-sequence FGK-stars is very rare (see e.g. 
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ICarpenter et alil2009l) . Since dust at similar radii would emit at 
shorter wavelengths for earlier-type stars, and since we were un- 
able to establish a statistical difference in emission at the same 
wavelength, we conclude that we cannot establish a statistical 
difference for dust at similar radii either We note, however, that 
the exact behavior of dust around M-dwarfs is complicated and 
potentially different from the evolution around earlier-type stars, 
and that a detailed analysis would require modeling the dust pop- 
ulation. 

Some detections of very hot dust very close to the star emit- 
ting in the band have been made thro ugh interferometry (e.g. 
lAbsil et al.L 120091; IdiFolcoet all l2007h . However, this dust is 
close to the subUmation limit and might represent a very differ- 
ent dust population than what we probed with our study. 

4.3. Mass limits 

The mass limit for the dust is coupled to the total fractional ex- 
cess (together with the albedo, this can be converted into the 
total optical depth around the star) as well as the distance from 
the star The closer the dust is to the star, the lower the total dust 
mass. 

As pointed out above, the exact location of the dust is dif- 
ficult to calculate. However, we note that in general, we are 
more sensitive to small a mounts of dust than comparable stud- 
ies for earlier-type stars. iGautier et alJ (l2007h pointed out that 
their 24 yum observations of comparable depth led to upper lim- 
its for the dust mass as low as lO^'Mg, although they did not 
state their assumptions. Because the location of blackbody dust 

1 /2 

grains scales as Rdust ~ , the dust mass detection limit scales 
as Mdust ~ Li,, so the main reason we are able to place very 
strict limits on the dust masses is the low luminosities of the 
M-dwarfs. 

Assuming blackbody grains of Bond albedo 0.3 and a sin- 
gle grain size of 1 fim with a density of 2.5 g/crn^, we arrive 
at upper limits for the dust mass {3cr) of 1.7 * lO^^Me for the 
0.0035 Lq Bai-nard's star and 2.5 * 10 ''Me for the 0.01 1 L© stai- 
GJ 526 for dust of ~ 150K, and even stricter limits for dust at 
higher temperatures. We cannot make statements about dust at 
low temperatures, large amounts of which could be hidden from 
our survey at larger orbital radii. However, we are not aware of 
any detections of far-lR or sub-mm excesses for these objects. 

5. Implications 

The dust evolution around M-dwarfs is still significantly less 
constrained than the evolution around earlier- type stars. The evo- 
lution of dust might follow a different path due to the lower 
luminosity combined with the low surface temperatures (and 
thus peak emission at longer wavelengths, longer than the small- 
est grain sizes). Blowout seems t o be irrelevant f or M- dwarfs. 
Theoretical models developed by iPlavchan et al.l (|2005) imply 
that the fractional infrared luminosity L/^/L* might be indepen- 
dent of age for old M-dwarfs and there might be a fractional 
excess of ~ 9 * 10"^ for all M-dwarfs. It is important to point 
out that we would not be able to detect such an excess. We fur- 
thermore realize that an equal excess for all stars in our sample 
would not be detected at all - because it would be masked by our 
fitting routine. 

More investigations will be required to understand M-dwarf 
debris disks in more detail. While the Spitzer 24fj.m channel is 
not available any longer and it will be quite a while until we can 
receive data from JWST, there are other promising techniques 



available to find debris disks around M-dwarfs. For example, 
Herschel has strong capabilities in the far-infrared, which al- 
lows one to study the frequency of colder dust. Studies at even 
longer wavelengths are possible (and have been conducted) in 
the sub-mm regime. To arrive at a complete picture of the M- 
dwarf debris disk population, both short- and long-wavelength 
studies will be required. 

There are other interesting techniques on the horizon as 
well. Up to now, infrared excess studies have always been 
more sensitive to circumstellar dust than resolved imaging and 
therefore the amount of debris disks detected by their SED 
excess emission is much larger than the number of resolved 
debris disks. With the upcoming next-generation VLT instru- 
ment SPH ERE and its high-con trast polarimetric imager called 
ZIMPOL, dSchmid et al.L l2006V this could change. Owing to 
the use of the pola rimetric differential imaging technique (e.g. 
lOuanz et al.L l201ll) . combined with an extreme adaptive optics 
system, ZIMPOL will have a very small inner working angle 
(IWA) of ~ 0.02". Because a dust ring this close to the star only 
occupies a small area on the sky, the surface brightness is high 
even when the total luminosity of the disk in scattered light is 
very low. 

Assuming an unfavorable face-on viewing angle, blackbody 
grains, a thin dust ring compared to detector resolution, a geo- 
metric albedo of 0.3 and scattering albedo * polarization fraction 
- 0.1, we estimated the performance of SPHERE in detecting 
such rings. We converted these predictions into detection limits 
in terms of F/s/F* and plot these limits in figure|5] We used pre- 
liminary performance estimates for ZIMPOL (achievable con- 
trast of 7.8 / 8.1 / 9.5 (mag/arcsec^ - mag) at separations of 
0.02 / 0.2 / 1 "). SPHERE is expected to be more sensitive to 
warm dust than the analysis we present in this paper. The steep 
cutoff toward higher temperatures arises from the inner work- 
ing angle of ZIMPOL. As we discussed in section 14.21 the real 
dust of a given temperature might be farther out than expected 
from a blackbody grain calculation because of the grain size dis- 
tribution. Therefore, this might be able to probe warmer dust. 
Furthermore, this technique should work better (up to higher 
temperatures) for nearby sources. For example, for Barnard's 
star at 1.83pc, we can probe warmer dust than for GJ 526 at 
a distance of 5.41pc. SPHERE will not replace surveys in the 
mid- and far-infrared to sub-mm, but allow complementary ob- 
servations of nearby sources, probing them for hot dust at small 
separations and attempting to resolve the dust spatially. It allows 
one to target individual stars to clarify uncertain detections or 
systems that are on the edge of having significant excess, such 
as the GJ 860 A/B system discussed in section [J.6.4l It will also 
be able to study known debris disks such as that of AU Mic or 
possibly AT Mic. 

Together, all these techniques will hopefully enable us to 
understand not only the M-dwarf debris phenomenon, but de- 
bris disks in general in a more detailed way, including their fre- 
quency, their evolution, and their structure. 



6. Conclusions 

We performed a volume-limited search for mid-infrared ex- 
cesses at 12 and 22//m using the RECONS database and the re- 
sults of the 2MASS and WISE all-sky surveys. Our main findings 
are: 

1 . It is possible to fit the infrared colors of M-dwarfs with a 
polynomial using either the absolute magnitude in the first 
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WISE band (3.4/im) or the V-Ks error as a proxy for stel- 
lar temperature. The resulting fit enabled us to search for 
excesses down to typically 14.9% for the [3.4]-[12] color, 
17.9% for the [3.4]-[22] color and 7.0% for the [12]-[22] 
color (3cr median values). 

2. The quality of the WISE data is excellent. The reduced 
goodness-of-fit values for the polynomial fits suggest that 
either the WISE errors are overestimated or that there is a 
significant portion from a systematic calibration uncertainty 
in these errors, which would be corrected for by our fitting 
routine. 

3. No excess from a debris disk is found in any color within our 
detection limits. One physical excess is detected around SCR 
1845-6357 A, but is shown to come from the T6 companion 
orbiting this late M8.5 dwarf. This proves the validity of our 
method. We also found a tentative excess around AU Mic at 
22yum, but it is not part of our statistical sample. 

4. Compared to surveys for earlier-type stars, our derived ex- 
cess fraction of 0;^,',- y% is lower, but this is probably an age 
efi'ect. The excess fraction in the mid-infrared (warm debris) 
is not significantly difl^erent from the excess fraction for old 
FGK dwarfs. 

5. Dust evolution around M-dwarfs is complicated and compar- 
isons to earlier-type stars at either the same orbital distance 
of the dust or the same total dust mass are difiicult because 
the primary clearing mechanisms for small dust are difl^erent. 
We note, however, that the limits on the total dust mass that 
can be derived from our survey are more stringent than the 
typical limits for earlier-type stars. 

6. The upcoming SPHERE/ZIMPOL instrument wiU search se- 
lected targets for reflected light from warm, close-in dust. 
This technique should be more sensitive for nearby bright 
stars than either WISE or Spitzer and could provide a mech- 
anism to probe dust at very smaU angular separations. 

We conclude that WISE is a great tool for studying the in- 
frared properties of M-dwarfs and stars in general and will in the 
future provide a great source of information, just like IRAS did 
almost three decades ago. The detection limits achievable with 
WISE are close to the limits achievable with pointed observa- 
tions from Spitzer. 

In general, there is still a severe lack of knowledge about 
dust around M-dwarfs and more investigation is required both 
theoretically and observationally to understand the similarities 
and dififerences of debris dust evolution between M-dwarfs and 
earlier- type stars. While instruments like Herschel can help us 
to learn more about dust at far-infrared wavelengths, WISE and 
SPHERE/ZIMPOL can help by probing M-dwarfs for closer-in, 
warm dust rings. 
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Fig. A.l. Obtained fits and error estimates for all three colors and 
additionally the Ks - [22] and Ks - [3.4] colors against absolute 
magnitude in the first WISE band and against V-Ks color 

Appendix A: Details on the fitting routine 

For the determination of a fit curve to the color-magnitude data 
required for the subsequent search for excess sources, we used a 
polynomial fitting routine. We chose a polynomial degree n and 
then proceeded as follows: 

1 . Select only the single stars from the sample of all stars 

2. Perform a least-squares fit for polynomial of degree n 

3. Calculate goodness-of-fit 

4. Calculate cr' - for all stars 

5. Remove stars with |cr'| > 3 from the sample 

6. Repeat from 2) until converged 

We did this for polynomials of increasing degree, starting 
with n = 1 . For all colors, the quality of the fit does not increase 
significantly beyond n = 4, therefore we chose to use polynomi- 
als of fourth degree in all cases. 

From the value of the fit, it is possible to estimate the 
amount of intrinsic scatter in the stellar colors, assuming the er- 
ror estimates for WISE are correct. As seen before, the WISE 
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Fig. A.l. Inferred errors for our fits. As can be seen, the fit is very 
well constrained for intermediate M-dwarfs, but less constrained 
at both ends of the fit. At the high-luminosity end, we mitigate 
this problem by including several K-dwarfs. 



errors are likely to be overestimated because the values for 
all three fits are lower than one. We conclude that the intrinsic 
scatter is small and the actual deviations from the fit are dom- 
inated by measurement errors. This is true for all WISE colors 
([3.4]-[12], [3.4]-[22] and [12]-[22]). The K, - [22], which we 
also investigated and show here for reference, has a higher 
value, and intrinsic scatter (induced by scatter in the Ks color) 
might play a more substantial role here. Ks - [3.4] is very good 
again, with a low x^ value. We show it here for reference as well 
because it bridges the 2MASS longest and WISE shortest wave- 
lengths for M-dwarfs, which we think could be useful. 

While it is possible to obtain uncertainties on the derived 
parameters of a polynomial fit, this is of limited helpfulness be- 
cause the individual parameters are highly correlated. To esti- 
mate the error on our fit, we again took a data-driven approach. 
We used the unmodified cr errors from WISE (not cr') and pro- 
duced 1000 data sets by randomly adding Gaussian errors with 
a standard deviation of cr, to each measurement. We fit these 
1000 data sets using our fitting routine. We then calculated the 
uncertainty of our fit at every point by calculating the standard 
deviation of these 1000 obtained polynomial fits. 

Figure I A. II summarizes the fits and their upper and lower 
limits (dashed lines). Because the errors are small and hard to 
see in this representation. Figure rA2l shows the error of each fit 
as a function of the absolute magnitude in the first WISE band. 
As can be seen, the errors are generally small, mostly below 0.02 
magnitudes, and thus much smaller than the measurement errors 
from the individual stars, which is why we did not provide an 
in-depth discussion of fitting errors in our analysis section (the 
error is not dominated by the fitting error). At both ends of the 
fit, the fitting errors naturally increase because the fact that the 
fit is less constrained at the ends. For the high-luminosity end 
of the M-dwarfs, this effect can be mitigated by including some 
K dwarfs, as we have done. We have good reasons to believe 
that our fundamental assumption that the colors can be fitted by 
a (fourth-order) polynomial holds reasonably well in the stellar 
regime (namely, the;^'^ values of our fits are very good). 

We summarize the parameters of our obtained fourth-order 
polynomials in table lA.ll The expected color can then be calcu- 
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Table A.l. Polynomial fits 



Color 


Proxy 










CI4 




41 


- ri2i 


M3.4 


5.6839e+00 


-3.9338e+00 


9.5780e-01 


-9.7494e-02 


3.5935e-03 


0.6943 


[3.4] 


-[22] 


M3.4 


8.2659e+00 


-5.6808e+00 


1.3899e-^00 


-1.4269e-01 


5.3063e-03 


0.6054 


[12] 


-[22] 


™3.4 


i . / JO lc~rUU 


-1.1/ Z Ici-UU 


Z.077IC-UI 




1 1 9QSfi 0^ 




Ks- 


-[22] 


M3.4 


6.6826e-^00 


-4.7923e-t-00 


1.2182e-f00 


-1.2829e-01 


4.8876e-03 


1.6584 


Ks- 


■ [3.4] 


M3.4 


2.2181e-f00 


-1.4508e-t-00 


3.4410e-01 


-3.4285e-02 


1.2482e-03 


0.5662 


[3.4] 


-[12] 


V-K, 


5.5283e-01 


-7.6010e-01 


2.9965e-01 


-4.1715e-02 


1.9851e-03 


1.0173 


[3.4] 


-[22] 


V-K, 


7.3936e-01 


-l.OlSle-t-OO 


4.2361e-01 


-6.2772e-02 


3.1802e-03 


0.7858 


[12] 


-[22] 


V-K, 


2.7830e-01 


-2.9003e-01 


1.1597e-01 


-1.7087e-02 


8.6622e-04 


0.5322 




-[22] 


V-K, 


1.261 le-f 00 


-1.491 le-l-00 


5.8420e-01 


-8.3413e-02 


4.1207e-03 


2.4038 




■ [3.4] 


V-K, 


1.2223e-f00 


-1.1065e-i-00 


3.5304e-01 


-4.4812e-02 


2.0075e-03 


0.4231 



Notes. Parameters of our estimated polynomial fits for the different colors and the different stellar temperature proxies (3.4yum absolute magnitude 
and V-ATj color). Ax~ goodness-of-fit value is given for reference. Description of the parameters in the text. 

lated as 

color(p) - flo + a\p + a2p^ + a^p^ + a^p'^ , 

with p being the proxy for the stellar temperature (either absolute 
3.4//m magnitude or V-Ks color). In addition to the three WISE 
bands, we also provide fits for the Ks - [22] color here. The Ks- 
[22] color can also be used for Spitzer data using the small offset 
of 0.083 magnitudes determined in Fig.|2l 
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